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ABSTRACT 

The nonlinear evolution of matter clustering induces temperature anisotropies in the cosmic microwave 
background. The fluctuations in temperature are produced at a relatively low redshift, z ^ 10, and con- 
sequently are not affected by the uncertain reionization proccesses that occurred at earlier epochs. The 
amplitude of the effect depends on the evolution model assumed for galaxy clustering. For plausible param- 
eter choices, the temperature fluctuations are in the range 10~^ ^ AT/T ^ 10^^. Observed limits on AT/T 
would be violated if nonlinear clustering on present cluster scales commenced prior to z ^ 10. 

Subject headings: cosmic background radiation - cosmology 
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I. INTRODUCTION 



Most of the studies on the anisotropy expected in the temperature of the Cosmic Microwave Background 

(CMB) have been based on linear density perturbations. Comparatively little attention has been given 
to additional effects that can be generated during the nonlinear phase of matter clustering evolution: the 
Sunyaev-Zeldovich effect relevant on small angular scales (Scaramella et al. 1993), bulk motions of gas (Os- 
trikcr and Vishniac 1986; Vishniac 1987; Hu, Scott and Silk 1994) and time-varying gravitational potentials. 
The formalism to calculate the temperature ansotropies generated by time-varying potentials during the 
nonlinear phase has been developed elsewhere (Martinez-Gonzalez et al. 1990) for a flat background. Pro- 
vided that one considers scales much smaller than the horizon and non-rclativistic matter peculiar velocities, 
it is possible to describe the gravitational field in terms of a single potential 0(i; x) that satisfies the Poisson 
equation 



where a{t) represents the scale factor, x are comoving coordinates, pi, (x a ^ is the background density and 
A is the density fluctuation defined by p = pb{l + A). We have chosen units such that SttG = c = 1. 

The temperature anisotropics are given by the expression 



where n is the unit vector in the direction of the observation and subscript LS denotes the last scattering 
surface. In linear theory, the gravitational potential due to the linear density fluctuations is static if = 1, 
so only the two first terms survive (Sachs and Wolfe 1967) to give the primary anistropies. Unless il < 1, the 
integrated effect along the path of the photon only contributes as the fluctuations enter the non- linear regime, 
to induce secondary anisotropics. The effect of second-order perturbations on the CMB has been previously 
studied (Martinez-Gonzalez et al. 1992), and were found to contribute to AT/T at a level of ^ 10^® for 
popular models of galaxy formation. The effect of isolated structures (voids, clusters, great attractors, etc.) 
has also been estimated to be at a level of AT/T - 10""^ - 10"*^ (R ecs and Sciama 1968; Thompson and 
Vishniac 1987; Kaiser 1982; Nottale 1984; Martinez-Gonzalez and Sanz 1990). Little attention has been 
paid to the intermediate regime of weak clustering. The effect of clustered matter on the CMB in the very 
non-linear regime has been estimated using N-body simulations, for isolated structures (van Kampen and 
Martinez-Gonzalez 1991) and for the hot dark matter (HDM) scenario (Anninos et al 1991), where effects of 
the order of 10^® and 10^'"' are found, respectively (in the latter paper, they also include Tompson scattering). 



We consider a flat Friedmann dust model to represent our universe from recombination (z ~ 10^) to the 
present time {z = 0). The microwave photons are propagated since the last scattering surface {z ~ 10^) 
in any direction and are aff'ected by the gravitational potential due to the density fluctuations which are 
evolving from the linear regime to the nonlinear one. We are interested in the effect of nonlinear fluctuations 
generating anisotropics on the CMB. We calculate this effect assuming a 2-point density correlation function 
that can be justified from observations at the present time and take the time evolution from numerical 
simulations. 

We wish to calculate the effect of nonlinear density fluctuations operating at late times (typically z ^ Q) 
on scales much smaller than the horizon. In this case, a non-static potential <p{t, x) suffices to represent the 
gravitational field of the matter and its efii'ect on the photons (Martinez-Gonzalez et al. 1990), whence 



V2(/. = -p6a2A(t;x), 



(1) 




(2) 



II. THE NONLINEAR GRAVITATIONAL EFFECT 



a) The effect of a time-varying potential 




(3) 



V'^(t){t,x) = 6a"^A(t,x), 



(4) 
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where A represents the nonUncar density fluctuations. The integral (3) is extcmded along the photon path: 
x(t,n) = A(a)n (A(a) = 1 - a^) gives the distance to the photon from the observer normalized to the 
present horizon distance. All units have been chosen as in the previous section and the scale factor and 
horizon distance at the present time are unity, i.e. Uo = 3to = 1- 

The temperature correlation function C{a) is obtained from the previous equation by averaging over all 
direction pairs separated by an angle a 

C{a} = 4 J\x^ j\x,.^^^C^{X,,X^-x), (5) 

x(Ai, A2, a) = |Aini — A2n2| = [A^ + X^ — 2Ai A2 cosa]^/^, (6) 

where rii • n2 = cosa. Thus, C{a) is given in terms of the correlation of the gravitational potential 
C^(Ai, A2;a;) at two different cosmic times Ai and A2. 



b) Models for the evolution of galaxy clustering 

Wc do not have information from observations about the matter correlations at different cosmic times. 
Observations of the galaxy and galaxy cluster distributions are limited to the local universe with approxi- 
mately the same cosmic time for all galaxies and clusters. However, numerical simulations of current models 
of galaxy formation (Cole and Efstathiou 1989; Hamilton et al. 1991) imply that there is a simple time 
scaling for the evolution of the matter correlation function. We will assume the following 2-point correlation 
function ^(^1, 2:2; r) 

^{zuZ2;r) = [{l + z,){l + Z2)]-'i'C-^)\ rKr^, 
= 0, r> Tm- 

Here r represents the physical coordinate and r„i is a large-scale cut-off beyond which any correlations are 
assumed to be vanishingly small. The previous formula can be justified as follows. The spatial dependence 
is obtained from present epoch observations of galaxy clustering, 7 ~ 1.8 and ro — 5/i~^Mpc (Davis and 
Peebles 1983; De Lapparcnt ct al. 1988), which extend to a scale r„i ~ 15/?,^^Mpc. However, positive 
correlations have recently been found in the APM catalog (Maddox et al. 1990) up to angular scales of ^ 5° 
(corresponding to length scales of 30/i^^Mpc, given the magnitude limit of the catalog). An equivalent 
signal is also apparent in the counts in cells of IRAS galaxies of amplitude AN/N 2± 0.5 in cubes of 30 /i^^Mpc 
size (Saunders et al. 1991). Additional evidence of clustering up to scales of ~ 50 h~^Mpc comes from the 
distribution of galaxy clusters (Bahcall and Soneira 1983, Dalton et al. 1992, Peacock and West 1992), 
although the amplitude of the cluster autocorrelations is strongly biased. 

The evolution of correlations in numerical simulations with scale-free spectra as initial conditions, P{k) oc 
fc"",— 1 < n < 0, can be approximately described by the power law £,{z,r) = (1 -|- z)^'^{ro/r)'^ for scales 
r < 2ro(l + z)~^/'^ (Cole and Efstathiou 1989) at late times, {z ^ 6). This is equivalent to setting e = in 
(7) and describes the evolution of stable clustering where clusters are already formed and virialized at high 
z. Less radical evolution of clustering occurs if C(r, z) is constant with time in comoving coordinates. In this 
case, the clustering developed early and froze out, as would be appropriate for a biased galaxy formation 
model. In this case,, e = 7 — 3. It is also possible to have e > in models where clusters form more recently 
and are still evolving at present. Some N-body simulations (Melott 1992) predict < e < 3 for scenarios 
with spectral indices in the range — 3 < n < 1. 

Wolfe (1993) has recently found that the clustering of damped Lya absortion systems (2 « 2 — 3) in 
quasars is best fitted with an evolution exponent e = —1.2 (7 = 1.8). In what follows, we will not choose 
any specific value of e, but we will consider values within a broad interval (—1.2 < e < 2) that covers the 
range of interest. Roughly speaking, values of7 — 3<e<0 will represent different evolutionary models 
where clusters form very early, whereas e > describes late and ongoing clustering, as would be appropriate, 
for example, if O << 1 and = 1, respectively. 
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The above expression for ^ can be rewritten in terms of A and comoving distance x as 



^(Ai, A2; x) = [(1 - Ai)(l - A2)]3-^+^(^r, X < Xm, 

= X>Xm- (8) 

We will test this type of 2-point correlation function in terms of the anisotropics generated in the CMB. 
The ad hoc cut-off Xm is time dependent and represents the maximum scale where the nonlinear evolution 
model is appropriate: ^(Ai, A2; x) — £^c', ~ 0.1. We write 

a:™ = ro$c"^[(l-Ai)(l-A2)]^. (9) 

Thus wc have three phenomenological parameters to describe clustering evolution: the initial rcdshift for 
the start of nonlinear evolution Zm (Am, = 1 — (1 + Ze)^^^^}, the evolution parameter e and the correlation 

cut-off Xm at present, Xmo = fm = fo(,c ^ • 

c) AT/T for models of galaxy clustering 
From the Poisson equation (4) one obtains 

V^, V^, C^(Ai, A2; x) = 36 [(1 - Ai)(l - A2)]-'e(Ai, A2; x) , (10) 

Assuming standard boundary conditions: i) jC^I < 00, ii) 6*0(00) — 0, iii) C^(0) = and iv) continuity 
up to the second derivative of C^, the correlation of the gravitational potential, associated to the correlation 
of matter given by equation (8), is unique and defined by the expression 

= /5^M-^{3[(4 - 7)(3 - 7) - 2] - (5 - - ^){ — f + 6( — )4-^}, x < x^, 

= /5x2a;i-''2(4-7)(2-7)(^), x > x^, (11) 

where 

/ = [(1 - Ai)(l - A2)]i-^-, g = --^ + + (12) 

0—7 4—7 0—7 2—7 

In the next section we illustrate these results for the potential correlation , C^{x), and the temperature 
correlation, C{a) (as given by equations 5 and 11), for several cases of galaxy clustering evolution. 

III. RESULTS 



The behaviour of the potential correlation function C0{x) is shown in Fig. 1 as a function of the comoving 
distance x for several values of the evolution parameter e and intial redshift Zm- The correlations grow with 
z for comoving clustering (e = —1.2) whereas they decrease for stable clustering (e = 0) and for more rapid 
evolution (e > 0). This behaviour of implies that the nonlinear effect on AT/T depends on the initial 
redshift at which the nonlinear evolution of matter initially occurs for comoving clustering. For stable and 
faster clustering evolution, the effect will be produced at low redshift z ^3 even if we extrapolate the power 
law evolution (cq. 8) to arbitrarily high z. 

Fig. 2 shows the temperature correlation function generated by the nonlinear evolution of matter given 
by eq. (8) for three cases of the evolution parameter e = —1.2, 0, 2. The correlation scale is of the order of a 
degree and is larger for higher e. In the case of comoving clustciring (fig. 2a), most of the effect is produced 
close to the initial redshift z ^ 10. However for stable clustering the temperature correlations are generated 
in the range 3 ^ z ^ 0.5, (fig. 2b) and for e = 2 in the range 0.5 ^ z ^ 0.1 (fig. 2c). 

We also find that C{a) is very sensitive to the matter correlation cut-off ^c- This is shown in fig. 2d,e,f 
for = 0.287(r™ = 10/i-iMpc),0.2,0.1. 
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The main results that we have obtained for the different models of galaxy clustering evolution can be 

summarized as follows: 

(a) The temperature correlations C{a)i ^ 10~^ for values of the evolution parameter —1.2 < e < 3, initial 
redshift -Zm > 3 and cut-off > 10 h~^Mpc at present (i.e. < 0.287). 

(b) The effect on AT/T reaches a maximum for comoving clustering (e = —1.2), decreases for stable clustering 
(e = 0) and slowly increases for e > 0. 

(c) More than 80% of the nonlinear effect is produced in the interval z ^ 3 for e ^ whereas for e = —1.2 
the effect occiirs near the initial redshift when eomoving nonlinear evolution commences. 

(d) Considering present experimental limits on degree scales, the evolution of the correlations for comoving 
clustering as given by eq. (8) (with = 10/i~^Mpc) must have commenced at 2; ^ 10, otherwise the 
present limits on AT/T of w 10~^ would be violated. Extrapolating the observed galaxy correlation 
power law up to Vm = 19.4/i~^Mpc (^c = 0.1) those limits imply that clusters must have formed later 
than 2: ~ 4 for the case of comoving evolution. 

IV. CONCLUSIONS 

We have shown that galaxy clustering evolution typically generates CMB temperature fluctuations in the 
range 10~^ ^ AT/T ^ 10~^. For some specific models, an observable effect on degree scales can be obtained. 
This anisotropy is generated at redshifts below ~ 10 and therefore is not erased by reionization processes 
that could take place in the universe at high z. 

Observational limits on the anisotropies of the CMB on degree scales constrain models for the evolution 
of galaxy clustering and therefore of cluster formation. The maximum anisotropy is generated in the case of 
comoving clustering where the limits on AT/T are violated if we extrapolate the observed galaxy correlation 
power law up to 19.4 h~^Mpc (corresponding to a correlation = 0.1) with Zm = Q or allowing Zm ^ 10 and 
keeping the cut-off at 10 h~^Mpc (^c = 0.287). The stable clustering model (e = 0) or models with a faster 
evolution (e > 0) contribute to AT/T at a level of the order of 10~^ whereas models that form structure at 
high-z are constrained by the CMB observations. 

It would be very interesting to calculate the nonlinear gravitational effect using more realistic models for 
the evolution of galaxy clustering that smoothly interpolate between the linear and nonlinear regimes. Such 
a calculation would clarify the nonlinear versus linear contributions. Indeed, from the numerical simulations 
done by Anninos et al. (1991) for a HDM model, it is seen that the nonlinear gravitational effect can be of 
the same order as the linear Sachs- Wolfe effect on degree angular scales. 
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FIGURE CAPTIONS 



Figure 1. Potential correlation function C^{x) as a function of the comoving distance x for three values 
of the evolution parameter: a) e = —1.2 (comoving clustering), b) e = (stable clustering) and c) e = 2. 

Different lines correspond to the following rcdshifts for the start of nonlinear evolution: = (solid), 0.5 
(dotted), 1 (short-dashed), 3 (long-dashed) and 6 (dash-dotted). The correlation cut-off is assumed to be 
= 0.287. 

Figure 2. Temperature correlation function C(a) as a function of the angular distance a. a) e = —1.2 
and Zm = 10,6,3,1 (solid, dotted, short-dashed, long-dashed), b) e = and Zm = 6,3,1 (solid, dotted, 
short-dashed) and c) e = 2 and Zm = 6.0.5.0.1 (solid, dotted, short-dashed). The correlation cut-off is 
assumed to be = 0.287 for cases a,b,c. d) e = —1.2, e) e = 0, f) e = 2, and = 0.287, 0.2, 0.1(rTO = 
10, 12.2, 19.4 /i~^Mpc) for solid, dotted and dashed line respectively. The starting redshift 2:^ = 6 is the 
same for cases d,e,f. 
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